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The parton momentum densities in the proton are determined from a NLO QCD analysis of structure functions 
measured by HERA and fixed target experiments. Also included are data on the difference of the up and down 
anti-quark densities. The uncertainties in the parton densities, structure functions and related cross sections are 
estimated from the experimental errors, taking into account all correlations. Standard Model predictions of the 
charged current Born cross sections at large x and Q 2 are calculated and compared with recent data from ZEUS. 



1. INTRODUCTION 

With the integrated luminosity of about 
50 pb _1 collected at HERA during the years 
1994-1997 a new kinematic domain of large x and 
Q 2 becomes accessible for the study of deep in- 
elastic scattering in ep collisions. Measurements 
by ZEUS of the e+p single differential neutral cur- 
rent (NC) and charged current (CC) Born cross 
sections for Q 2 > 200 GeV 2 have recently become 
available [@,§. 

Standard Model (SM) predictions calculated 
with, for instance, the parton distribution set 
CTEQ4 [ H are in good agreement with the NC 
cross sections but fall below the CC measure- 
ments at large x and Q 2 . This could be an indi- 
cation of new physics beyond the SM but might 
also be due to an imperfect knowledge of the par- 
ton densities in this kinematic region. For in- 
stance in [ it is shown that a modification of 
the CTEQ4 down quark density yields SM pre- 
dictions in agreement with the CC e + p data. 

To investigate these issues we have performed a 
global NLO QCD analysis [ || of structure func- 
tion data to obtain the parton densities in the 
proton. A full error analysis provides estimates 
of the uncertainties in these parton densities, the 
structure functions and the SM predictions of the 
NC and CC cross sections. 

2. QCD ANALYSIS 

The data used in the fit are F% from ZEUS [Jl| 
and HI [ § together with Ff , F 2 d [ |, § ffi| JL1 
and F2/F2 [ |l2| from fixed target experiments. 



Also included were neutrino data on xF^ Fe [ 13 
for x > 0.1 and data on the difference x(d — u) [ 
0- 

After cuts the data cover a kinematic range of 
10" 3 < x < 0.75, 3 < Q 2 < 5000 GeV 2 and 
W 2 > 7 GeV 2 . Both F 2 d and F$/F% as well as 
x pvFe were corrected for nuclear effects using the 
parameterization of [ [l5| for deuterium (with an 
assumed uncertainty of 100%) and that of [ [Tg| 
for iron (assumed uncertainty 50%). 

The QCD predictions for the structure func- 
tions were obtained by solving the QCD evolu- 
tion equations in NLO in the MS scheme [ [[7). 
At the scale Qq = 4 GeV 2 the gluon density (xg) , 
the sea and valence quark densities (xS, xu v ,xd v ) 
and the difference x(d — u) were parameterized in 
the standard way [ [l8| (16 free parameters). The 
strange quark density was taken to be a fraction 
K s = 0.20(±0.03) of the sea [ [jj]]. The normal- 
izations of the parton densities were constrained 
such that the momentum sum rule and the va- 
lence quark counting rules are satisfied. The 
charm and bottom quarks were assumed to be 
massless and were generated dynamically above 
the thresholds Q\ = 4(±1) and Q\ = 30 GeV 2 
respectively. The value of the strong coupling 
constant was set to a s (M§) = 0.118(±0.005). 

Higher twist contributions to F 2 and F$ were 
taken into account phcnomcnologically by de- 
scribing these structure functions as 

F »T = F LT [l + H{x)m (1) 

where F^ T obeys the NLO QCD evolution equa- 
tions and where H(x) was parameterized as a 
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fourth degree polynomial in x (5 free parameters). 

The normalizations of the ZEUS, HI and NMC 
data were kept fixed to unity whereas those of the 
remaining data sets were allowed to float within 
the quoted normalization errors (7 parameters) . 

The uncertainties in the parton densities, struc- 
ture functions and related cross sections were es- 
timated taking into account all correlations. In- 
cluded in the error calculation are the experi- 
mental statistical errors, 57 independent sources 
of systematic uncertainty (propagated using the 
technique described in [pTJ) and the errors on the 
input parameters of the fit. 



Eq. ([I]) is found to be negative for x < 0.5 and 
becomes large and positive at high x. Our result 
on the higher twist contribution is very close to 
that obtained by MRST [||. 

Nuclear effects in neutrino scattering were in- 
vestigated by calculating the ratios xF% Fe /xF^ N 
and FZ Fe /F£ N of the CCFR data to the QCD 
predictions of scattering on a free nucleon. Both 
these ratios do not significantly depend on Q 2 in 
accordance with measurements of nuclear effects 
in muon scattering [ ^2|. The ratios, averaged 
over Q 2 , are plotted in Fig. || and show the typi- 



3. RESULTS 



The fit yielded a good description of the data 
with a x 2 °f 1540 for 1578 data points and 28 free 
parameters. The quality of the fit is illustrated in 
Fig. [j] where we show the fixed target F% data 
for x > 0.1. The higher twist coefficient H(x) in 
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Figure 1. 



Fixed target data on F% versus Q 2 for 



x > 0.1. The full (dotted) curves correspond to the QCD 
prediction with (without) higher twist contributions. 
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Figure 2. Nuclear effects in u-Fe scattering estimated 
from a comparison of the CCFR data and the QCD fit 
as explained in the text. The curves correspond to the 
parameterizations from [ [lfij (used to correct xF^ Fe for 
x > 0.1). The errors from the QCD fit are drawn as bands 
around the curves. 



cal x dependence of nuclear effects, including the 
rise at large x due to Fermi motion. However, 
within the present accuracy it is not possible to 
establish whether these nuclear effects are signif- 
icantly different from those measured in charged 
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lepton scattering (full curves in Fig. |J). 

In Fig. H are shown the parton densities ob- 



They find B = 0.10 ± 0.01 which implies that 
d/u — > 0.2 as x — y 1. The d/it ratio is shown in 
Fig. U It is seen that the result from the QCD 



X9i 



: QCD fit 


Q 2 = 1 GeV 2 


- CTEQ4M 




& CTEQ5M 




s \ xg (x o.05) 


/ \ 


% " -0\ A 


3 xd v \ 


^xS (X 0.05) y^ffl^ 







io -3 icf 2 to - ' X 



d/u 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 X 



Figure 3. The parton densities from this analysis 

compared to the results from CTEQ4 and CTEQ5. The 
hatched bands indicate the error on the QCD fit. 



tained from the QCD fit (full curves). We have 
verified that the strange quark density is compat- 
ible with the measurements from CCFR [ |l~9| and 
also that the QCD prediction of F% agrees well 
with the ZEUS data [ || above Q 2 = 10 GeV 2 . 
This supports the assumption made in this analy- 
sis that, at least for x > 10 -3 , quark mass effects 
do not spoil the QCD extrapolations to large Q 2 . 

Also shown in Fig. || are the parton densities 
from CTEQ4 and CTEQ5 [ There is, within 
errors, good agreement between these results and 
the QCD fit, although both the present analysis 
and CTEQ5 yield a slightly harder xd v density 
than CTEQ4. Bodek and Yang [ [| also obtain 
a harder xd density by modifying CTEQ4 (where 
the ratio d/u — > as x — * 1) such that 

d/u -> d'/u = d/u + Bx(l + x). (2) 



Figure 4. The ratio (d + d)/(u + u) from the QCD fit 
compared to results from CTEQ4, CTEQ5 and CTEQ4, 
modified according to Eq. 13) with B = 0.1. The hatched 
band shows the error on the QCD fit. 



fit with B = (full curve) is for x < 0.75 close to 
the modified CTEQ4 distribution with B = 0.1 
(dotted curve) ; if we leave B a free parameter in 
the fit, we obtain B = —0.02 ± 0.01 (statistical 
error), close to zero. In any case, the large error 
band clearly indicates that the exact behavior of 
d/u at large x is not well constrained. It might go 
to zero (CTEQ), to a constant (Bodek and Yang) 
or may even diverge (this analysis) as x — > 1. 

The CC e + p — > uX cross section is predomi- 
nantly sensitive to the d quark density 

d 2 a / dxdQ 2 oc (1 — y) 2 (xd + xs) + xu + xc (3) 

where y = Q 2 /(xs) with s the ep center of mass 
energy (-10 5 GeV 2 at HERA). In Fig. | we show 
the ZEUS measurements [ ^ of the CC e + p single 
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Figure 5. The CC e+p cross sections dcr/dQ 2 (top) 
and dcr/dx (bottom) measured by ZEUS, normalized to 
the QCD predictions based on CTEQ4. The full (dotted) 
curves are the predictions from the QCD fit (CTEQ5). 
The hatched bands show the error on the QCD fit predic- 
tion. 



differential cross sections da/dQ 2 and da/dx for 
Q 2 > 200 GeV 2 , normalized to the NLO predic- 
tions calculated with the CTEQ4 parton distribu- 
tion set. The full (dotted) curves correspond to 
the predictions from the QCD fit (CTEQ5) which 
both achieve a better description of the data at 
large Q 2 and x due to an improved determina- 
tion of the d quark density. SM predictions of the 
e + p NC cross sections, calculated with the parton 
densities obtained in this analysis, also agree very 
well with the recent ZEUS data [ |l| (not shown) . 

We conclude that, within the present exper- 
imental accuracy, no significant deviations can 
be observed between the data and the Standard 
Model predictions. 
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